Real-time phasefront detector for heterodyne interferometers 
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We present a real-time differential phasefront detector sensitive to better than 3 mrad rms, which 
corresponds to a precision of about 500 pm. This detector performs a spatially resolving measure- 
ment of the phasefront of a heterodyne interferometer, with heterodyne frequencies up to approx- 
imately 10 kHz. This instrument was developed as part of the research for the LISA Technology 
Package (LTP) interferometer, and will assist in the manufacture of its flight model. Due to the 
advantages this instrument offers, it also has general applications in optical metrology. 
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I. INTRODUCTION 

Many applications in optical metrology require pre- 
cision measurements and characterization of laser beam 
wavefronts, as well as an accurate mode-matching of laser 
beams. To this end, it is usual to perform various adjust- 
ments that are both complex and time-consuming, based 
on repeated measurements of beam parameters [T] . Al- 
ternative methods are Shack-Hartmann sensors, which 
measure the shape of a single wavefront with an accu- 
racy of typically A/100 [2], and phase-shifting interfer- 
ometry [3JI1] (PSI), which reaches typically A/ 1000 [5]. 
Similar to PSI, the relative phase is found by a math- 
ematical algorithm [51 [7] which is applied on the inten- 
sities sampled for n > 4 instantaneous operating points 
with equidistant phase increments. This happens in par- 
allel at each pixel of the spatially resolving photodetec- 
tor (CCD camera). The advantage of our method over 
conventional PSI is that the phase-shift, normally imple- 
mented by a piezo-electric transducer (PZT), is intrinsi- 
cally contained in the time-dependent sinusoid obtained 
from the interference of two electric fields at different 
frequencies (optical heterodyning). Thus, non-linearities 
of PZT-elements, and additional electronics required to 
keep the interferometer at a specific operation point can 
be avoided [5] . 

This article describes an apparatus able to measure in 
real-time the relative spatial structure between two in- 
terfering wavefronts in a heterodyne interferometer. This 
allows a quick and very precise characterization of the 
mode-mismatch between two beams in real-time, making 
possible an online adjustment of the optical components 
according to the mode-mismatch displayed. By using 
beams with known shape, this method can also be used 
to analyze surfaces and optical components (in transmis- 
sion and reflection), reaching a sensitivity better than 
A/2000, which corresponds in our case to a precision of 
about 500pm at A = 1064 nm (Nd:YAG laser). 
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Optical heterodyne interferometry [9l \10\ is a useful tech- 
nique to measure distance variations with sub- wavelength 
precision and large dynamic range. This concept is 
applied, for example, in the LISA Technology Pack- 
age (LTP) QT] which utilizes a set of heterodyne Mach- 
Zchndcr interferometers to measure relative changes in 
the separation of two drag-free test masses with a noise 
level better than 10pm/vHz in the frequency range of 
3 mHz to 30 mHz. It is well-known that heterodyne inter- 
ferometers are susceptible to various noise sources [T21 03] 
which corrupt the phase measurement. One important 
noise source is the wavefront imperfections of the inter- 
fering beams in combination with beam jitter [T4l 115) . 
This effect is of particular importance if quadrant photo- 
diodes (QPD) are used, which are often employed to ob- 
tain alignment signals from the interferometer. The error 
term induced by the spatial inhomogeneity of the wave- 
fronts can be minimized if the interfering beams have 
identical shape. The device we present here allows one 
to match the shape of the wavefronts with a simple pro- 
cedure, permitting an online adjustment of the optical 
components. 



II. THEORETICAL BACKGROUND 

A. Interference pattern and heterodyne 
interferometry 

The electric field Ej (r,£) of a linearly polarized light 
beam can be described as 

E,- (r, t) = Ej p, exp { i [ 2tt /, t + <pj + i> s (r) ] } , (1) 



where j is an index to distinguish several beams, the vec- 
tor p ■ describes the polarization, Ej is the amplitude of 
the electric field, and ipj (r) is the spatial distribution of 
the phasefront. The intensity distribution I (r,i) oscil- 
lates at the heterodyne frequency /hot = f2 ~ fit an d 
is proportional to | E to tai(r,£) | 2 , where E to tai(r,£) = 
Ei (r, t) + E 2 (r, t) is the interference pattern. Assuming 
identical polarization vectors p x = p 2 and also that the 



differential coherence length of the laser used is much 
larger than the distance from the source to the recombi- 
nation point of the interferometer, the heterodyne com- 
ponent can be described as 

/ (r, t) = A{v) { 1 + C(r) cos [ 2-k f hct t - ip(r) ] } , (2) 

where A(t) is a space-dependent factor, C(r) describes 
the contrast in terms of the space r, and 



ip(r) = [ (pi + i>i (r) ] - [ <p 2 + if>2 (r) ; 



(3) 



gives the spatial dependence of the phase, where ip (r) = 
ipi (r) — ■02 (r) is the spatially distortion of the phase 
distribution. Ideally, for identically shaped wavefronts 
(01 (r) = -02 ( r ))i the longitudinal phase term (pi — ip2, 
which is the usual interferometric quantity to be mea- 
sured, contains the length measurement describing the 
pathlength difference AL between the arms of the inter- 
ferometer: 

where A = c/f is the wavelength of the light. This path- 
length difference can be strongly influenced by environ- 
mental changes that usually (as in LTP) disturb the main 
length measurement. 

In order to maximize the contrast of the interferometer, 
the beams need to be matched. It is therefore very use- 
ful to characterize the mismatch between them, and to 
be able to optimize it in real-time. 



B. Spatially resolved phase measurement 

The relative geometry of two interfering beams can be 
described by the spatial structure of the functions A(t), 
C(r), and ip (r). The apparatus described here measures 
these three functions in the real-time interference pat- 
tern using a CCD camera and pixelwise data processing. 
In order to obtain the phase of a sinusoidal function as 
in Equation pi) , several mathematical approaches [3] can 
be used which are based on measuring n equidistant in- 
tensity samples Ik — I{tk), with tk = kAt, where k 
is an integer. The approach chosen for this experiment 
is a 4-point algorithm |16j , corresponding to a straight- 
forward Discrete Fourier Transform (DFT) of the signal 
with n = 4 samples and At = T/4, where T = 1// is 
the period of the signal. For a noise-free signal, these 
intensities would be given by 



1 + C cos 



.(*+*!) 



(5) 



where J aV g is the average of the sampled intensities. 
The phase at the pixel 7, ip~, can be calculated from 



these intensity samples ( Iq 
algorithm as 

fry — arctan 



(7) 
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with the 4-point 
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Several other useful quantities can also be obtained from 
these 4 data points. Using the abbreviations 



a _j(7)_r(7) 

"7 ~ -'o 1 2 

°7 " 2 3 1 1 

d, = 1^+1^+1^+1^, 



(7) 
(8) 
(9) 



we get: 

• Contrast at the pixel 7: 



J<r:+^ 
° 7 - 2 ' d, 



• Total phase over the CCD surface: 

'/'total = arctan = 

V^7 a V 

• Total contrast over the CCD surface: 



(10) 



(11) 
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• Average intensity at the pixel 7: 

r(7) = ^1 
av e 4 



(12) 



(13) 



• The maximum and minimum intensity for the set 
of exposures (Iq . . . I 3) can also be determined for 
diagnostic purposes. 

Furthermore, an exposure of the dark fringe can be 
directly captured by triggering the CCD camera with 
the appropriate delay t^: 



37r/2 - (jCtotal „ 
Tdf — ~ 1 het • 



(14) 



III. 



INSTRUMENT COMPONENTS AND 
SETUP 



The interferometer configuration used in this experi- 
ment is a non-polarizing heterodyne Mach-Zehnder in- 
terferometer. The light source was a Nd:YAG NPRO 
(non-planar ring oscillator) laser with a wavelength of 
1064 nm. Two acousto-optic modulators (AOM), driven 
by slightly frequency shifted RF signals near 80 MHz, are 
used to generate two laser beams with a frequency differ- 
ence of /hot ~ 1623 Hz. Since the beams diffracted by an 
AOM have a distorted non-Gaussian beam profile, single- 
mode polarization-maintaining fiber optics are used as 
mode cleaners. The requirements on the CCD camera 
are: 



ij)(^ 




FIG. 1. Experimental setup used for the phasemeter. 

• Simultaneous exposure for every pixel ("global 
shutter" ). 

• The exposure time must be very short compared 
with the heterodyne period, Th t ~ 1/1623 Hz = 
616 fis. 

• The camera must be able to be triggered externally 
to allow frames to be captured at the required in- 
stances of time. 

• The signal for each pixel needs to be proportional 
to the intensity on that pixel. Saturation effects, 
such as "blooming" , must be avoided. 

The CCD camera used is a model XEVA-USB from Xen- 
ICs [T7j with a 12-bit dynamic range which reaches 30 
frames per second (fps) at a resolution of 320 x 256 pix- 
els with a pixel pitch of 30 fim. The photosensitive chip 
is made of InGaAs, which has a quantum efficiency of ap- 
proximately 80% for the near infrared (0.9-1.7/zm). The 
exposure time used was 1 fis (1/616 Xh c t). Ideally, the in- 
terference pattern should be sampled 4 times within a 
single heterodyne period of Th et . This would require a 
sampling period of At = 154 /zs (approximately 6500 fps) 
for the CCD camera which cannot be reached in practice 
due to the time required to transfer the image. In order 
to capture the intensity sample Ik, an integer number m 
of heterodyne periods Xh c t is added to At, and the trigger 
signal is sent to the CCD camera with the delay 



Atu 



mT] 



hct 



Tfc- 



(15) 



with Tfc = fcThot/4. The heterodyne period Xh e t is esti- 
mated by the timing control electronics of a Field Pro- 
grammable Gate Array (FPGA) at the beginning of the 
measurement as an average over 2000 periods, and is then 
transfered to the PC through the parallel port interface. 
The experimental setup is outlined in Figure [Tl In prac- 
tice, the environment is not stable enough to preserve a 
constant phase relationship over many periods of /hot- 
Hence, additional circuitry is used to re-synchronize the 
trigger timing electronics with the actual phase of the 
heterodyne signal: A single-element photodiode (SED) 
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FIG. 2. Time diagram of the signals processed to trigger the 
CCD camera. Note that the sinusoid is measured by the SED 
and yields a 180° phase shift wrt. to the interference pattern 
measured by the CCD camera. 



is located at the second output of the interference beam- 
splitter where the same interference pattern emerges with 
a 180° phase shift. The heterodyne signal measured by 
the SED is bandpass filtered and digitized by a compara- 
tor. When the phasemeter software is ready to capture 
the frame k, a command is sent from the PC to the FPGA 
through the parallel port, which includes the correspond- 
ing delay r^. The FPGA detects the rising edge of the 
digital heterodyne signal and triggers the camera with 
the controlled delay Tk (see Figure [2]). 

The CCD camera captures the frame and transfers 
it to the PC through its USB 2.0 interface. After all 
four frames have been acquired, the phasemeter soft- 
ware computes the physical quantities described in Equa- 
tions (I6J) to (13), and an additional exposure is captured 
by triggering the camera with a delay r<if given by Equa- 
tion (14), which corresponds to an exposure of the dark 



fringe. 



IV. RESULTS 



A. Real-time operation 



A graphical user interface (GUI) was developed to dis- 
play the measured data in real-time (see Figure |3| . The 
phasemeter reaches a rate of approximately 5 to 6 data 
displays per second, and hence enables real-time mode- 
matching of the beams, as well as an online optical align- 
ment of an interferometer and optimal adjustment of its 
components. The five different displays in Figure [3]are 
shown separately as 3D representations in Figures T|4] to 
[8] These measurements were conducted on a table-top 
Mach-Zehnder interferometer, as shown in Figure [I] 




FIG. 3. Graphical User Interface programmed to display the 
measured data in real-time. 
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FIG. 6. Exposure of a dark fringe. 
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FIG. 4. Spatial distribution of the phase. 
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FIG. 5. Spatial distribution of the contrast. 



FIG. 7. Average intensity over four exposures. 
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FIG. 8. Maximum intensity over four exposures. 



B. Performance of the phasemeter 

In order to measure the noise level of the instrument, 
the camera was illuminated with a spatially homogeneous 
light source (an array of infrared LED's behind a matt 
glass window), the intensity of which was modulated 
sinusoidally at a frequency of approximately 1623 Hz. 
This is a very stable test subject, ideal to determine 
the phase noise level of the instrument. It is expected 
to obtain a flat phasefront of the amplitude modulation 
from this measurement, since all pixels capture a simi- 
lar sinusoidally modulated intensity with the same phase 
relationship. A series of real-time phasefront measure- 
ments were conducted in this configuration, determin- 
ing the rms phase value at every pixel. Thus, a spatial 
rms phase variation of 2.96 mrad was meaured over the 
CCD area of observation. According to Equation Q, 
this value corresponds to a rms spatial resolution AL of 
about 500 pm. 

As it can be seen in Equation (pf, one noise source of 
the phase measurement is the fluctuation of the sampled 
intensities. The rms error of the phase, Aip Tms , induced 
by intensity fluctuations AJ rms can be estimated from 
Equation (pi) as: 



Ai^rms = , 2_^ 






A/ r , 



(16) 



(h - hf + (h - hf 



A/ r , 



After simplifying Equation (17) by using Equation ([5]), 
we obtain: 



Alarms = 



y/2 A/ rl 



C I, 



(17) 



ilVji, 



The following three error sources were identified, and 
their noise contribution to the phase measurement was 
estimated: 

1. Laser power fluctuations: An Allan deviation of 
8.6 x 10~ 4 was measured at an averaging time of 
33 ms, which corresponds to the sampling period of 
the CCD camera (30fps), yielding a phase error of 
1.22 mrad from Equation (Tl7|). 



2. ADC digital noise of the camera: A rms intensity 
error of four quantization units was measured by 
constant and spatially homogeneous illumination 
of the CCD camera. The pixelwise rms variation 
and an average over the CCD surface were then 
computed. This value corresponds to a relative 
intensity fluctuation of 9.76 x 10 -4 , which trans- 
lates (by using Equation (17) ) into a phase error 
of 1.38 mrad. 

3. Time jitter: There are at least three sources of jit- 
ter. Firstly, the synchronization delay of the com- 
parator output with respect to the 10 MHz clock of 



the FPGA, which is uniformly distributed between 
and 100 ns. Secondly, a similar delay between the 
FPGA clock and the CCD internal clock, which is 
also at 10 MHz but unsynchronized. Thirdly, other 
jitter effects such as apparent period fluctuations 
of the signal, due to setup conditional and limited 
stability, and the non-sychronization between the 
FPGA and CCD clocks with the clock of the mod- 
ulation electronics controlling the AOM's driving 
signal, and therefore the heterodyne frequency gen- 
eration. This latter effect is reduced (but not to- 
tally cancelled) by re-tracking the timing control 
electronics to the rising edge of the digital hetero- 
dyne signal for each exposure. A phase error of 
0.99 mrad was obtained by simulating the first and 
second effects in software, using two independent 
random delays uniformly distributed. 

Table U summarizes the noise contributions identified for 
the phase measurement. 



TABLE I. Main noise sources of the phase measure- 
ment. 



Noise Source 



RMS Phase Error 



Laser power fluctuations 1.22 mrad 

ADC digital noise of the camera 1.38 mrad 

Time jitter 0.99 mrad 

Total contribution 2.09 mrad 

Noise level measured 2.96 mrad 
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FIG. 9. Experimental setup with an additional lens in the 
path of one beam to intentionally change the curvature of its 
wavefront. 
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FIG. 10. (a) Phase front measured with a lens / = + 500 mm in one arm of the interferometer. The phasefront is clearly 
wrapped, due to the high curvature of the wavefront being transmitted through the lens wrt. the other one. (b) Phase front 
obtained by post-processing the data measured in Figure 10 [a) with a two-dimensional phase unwrapping algorithm. 
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FIG. 11. (a) Phase front measured with a lens f = — 500 mm in one arm of the interferometer, (b) Phase front obtained from 
post-processing the data of Figure [TT^a) with a 2D phase unwrapping algorithm. 



Measurements 



In order to test the functionality of the phasemeter, a 
series of measurements were conducted at the table-top 
interferometer shown in Figure IT] The aim of these mea- 
surements was to intentionally change the curvature of 



one of the interfering wavcfronts and to use the phaseme- 
ter to read out the spatial distribution of the phase. Since 
the two interfering wavefronts are similar, a lens was in- 
troduced into the path of one beam, between the corre- 
sponding fiber injector and the beam combiner (see Fig- 
ure I9J. The results of these measurements are shown 
in Figures [TU] to [12) Three different type of lenses were 
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FIG. 12. (a) Phase front measured with a cylindrical lens / = +80 mm in one arm of the interferometer, (b) Phase front 
obtained from post-processing the data of Figure 12 a) with a 2D phase unwrapping algorithm. 



used: / = +500 mm (Figure 10 1, / = —500mm (Fig- 
ure 11), and a cylindrical lens with / = +80 mm (Fig- 
ure 12). The curvature of the wavefront considerably 



increases (due to the lens) such that this covered sev- 
eral wavelengths in the area of observation. The result- 
ing phasefront measured is wrapped as can be recognized 
by the phase rings in Figures 10 'a) and 11 a), and the 
stripes in Figure [12^ a) . A two-dimensional phase un- 
wrapping algorithm [18tJ2TJ] was developed and was used 
for post-processing these data. The result of this post- 
processing is shown in Figures [lO^b ) , [TT^ b) , a nd[L2| ^b). It 
can be seen, by comparing Figures [J0[b) and|ll[bj, that 
the inflection of the phasefront curvature changes accord- 
ing to the focal length of the lens used (± 500 mm) re- 
spectively. Within the LTP interferometry research and 
development, this instrument will be used during manu- 
facture of the optical bench interferometer flight model. 
The lenses of the fiber injectors can be adjusted such that 
the difference in curvature between the two interfering 
wavefronts are minimized by using the real-time phase- 
front read out provided by this phasemeter. A phasefront 
measurement was already conducted at the optical bench 
engineering model for LTP [TTJ [ST] . These results are pre- 
sented in Figure 13 and clearly show an inhomogeneous 
phasefront, which can be attributed to non-optimal ad- 
justment of the lenses in the two fiber injectors. A fur- 
ther test was done on a table-top interferometer in order 
to adjust the lenses of two commercial fiber injectors by 
using this instrument. The aim of this adjustment was 
to match the parameters of the interfering beams and 
to obtain a homogeneous fiat phasefront. The result of 
this experiment is shown in Figure 14 A considerable 
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FIG. 13. Phase front measured at the engineering model of 
the optical bench for LTP. 



two beams (compare the phase scale of the plots in Fig- 
ures [l3] and IlP) . It can be seen that over a surface 



improvement was achieved in matching the shape of the 



of approximately 1.2 mm X 1.2 mm the phasefront shows 
a considerably homogeneous spatial profile. We analyzed 
a circular section of approximately 1 mm diameter at the 
beam center. The standard deviation of these data is 
3.49 mrad, which is very close to the measured sensitiv- 
ity of the instrument. This value corresponds to a spatial 
resolution of 590 pm, and amounts to a considerable im- 
provement in the correction of the beam shapes. 
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FIG. 14. Adjusted phasefront measured on a table-top Mach- 
Zehnder interferometer. 



V. CONCLUSIONS 

We have developed an instrument which allows real- 
time phasefront detection and mode-mismatch charac- 
terization of two interfering beams in a heterodyne in- 
terferometer. A rms noise level of 2.96 mrad, which cor- 
responds to a wavefront roughness of 500 pm, was ob- 
tained experimentally. This makes it possible to optimize 
the beam shapes by adjusting the optical components in 
real-time with the help of the data displayed onto the 
graphical user interface (5 to 6 data displays per sec- 
ond). The results shown in Figures 10 to 12 demonstrate 
the proper functionality of the instrument. By using well 
matched wavefronts, this method can also be used to an- 
alyze, and measure accurately, surfaces and optical com- 
ponents down to subnanometer levels. It is planned that 
this instrument be used in the manufacturing of the flight 
model of the optical bench for the LISA Technology Pack- 
age. 
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